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The charm squark resonance in the e^q channel provides one of the plausible interpretations 
of the reported anomaly at HERA. We show that the relevant R violating coupling A'^^gi 
is required to be large, typically around 0.1 in a large class of supergravity based models 
including the minimal one with the universal boundary condition at the GUT scale. Existing 
constraints on these couplings are reanalysed in this light a nd it is argued that such large 
couplings may be feasible but would require fine tuned cancelations. 



The anomalous events seen by the ZEUS and the HI detectors [^ at HERA in the deep 
inelastic e~^p scattering have generated considerable excitement [0-|^. These events would 
constitute evidence for physics beyond standard electroweak model if they are established 
firmly in the future. The presently available information when taken seriously allows two 
possible interpretations: (^) The presence of some lepton number violating contact interac- 
tion or (m) production of a resonance in the e^g channel [|]-P]. Supersymmetry, with 
violation of R parity |0-§1 provides a natural theoretical framework to incorporate the second 
possibility although an alternative in terms of a scalar leptoquark P] is also open. 

The supersymmetric interpretation of the HERA events assumes that the excess events 
seen at HERA are due to resonant production and subsequent decay of the squark to 
e~^q. Three possibilities have been considered in this context 0-§l: e^dji -^ ci^e^dji — >■ 
^L, e^Sij -^ ti- In analyzing these scenarios |0-§| it has been implicitly assumed that the 
squark masses are free parameters of the model. While this would be true in the most 
general situation, specific model dependence can alter some of the conclusions. Our aim 
is to show that the very minimal model dependent assumption on the charm squark mass 
necessarily requires large A'^21 ^o understand HERA events and this large coupling by itself 
is ruled out from other constraints. 

The specific assumption that we make and which leads to the above conclusion is that 
the charm squark mass squared is positive at the unification scale. This assumption is true 
in the radiative electro-weak breaking scenario with universal boundary conditions at the 
GUT scale, but it can also be true in a much more general context. We shall first assume 
that the gaugino masses are unified at Mqut but demonstrate later that the removal of this 
assumption does not significantly change the basic conclusion. 

The argument leading to the above conclusion is largely insensitive to the details of the 
radiative SU{2) x f/(l) breaking in the MSSM and runs as follows. 

Consider the following R violating couplings: 

^R = Kjki-^idiKij + eiUj)dl (1) 



The above terms are defined in the quark mass basis and K denotes the Kobayashi-Maskawa 
matrix. The charm squark interpretation of the HERA anomaly requires A'^21 to be non-zero. 
The HERA data can be explained provided 

0.025 - 0.034 

The number in the numerator of eq.(2) is indicative of the required range and depends upon 
the weightage given to the different experiments as well as on the next to leading order 
QCD corrections 0. In the following, we shall take the value 0.025 for the numerator 
in the RHS of eq.(2). B refers to the branching ratio for the squark decay to qe^ . This 
decay would take place through the coupling in eq.(l) itself. B is also influenced by the R 
conserving decays of the charm squark to an s (c) and a chargino (neutralino) . The A']^2i 
and the parameters /x, M2, tan/? determine B in the MSSM. HERA data can be reconciled 
if for a region in these parameters (z) eq.(2) is satisfied, (m) A'^21 i^ consistent with other 
constraints p!0|^l3| due to R breaking and [Hi) charm squark has a mass around 180- 220 



GeV. 

In supergravity based models, the charm squark mass at the weak scale is governed by 
the gauge couplings and the gaugino masses. Its value at Qo = 200 GeV is given in the limit 
of neglecting Kobayashi-Maskawa and cl — cr mixing by |14 

^L (Qo) ~ nil^ (Mgut) + 8.83M| + 1/2 M| cos 2(3 (1 - 4/3 sin^ Ow) (3) 

where we have assumed that the gauge couplings and the gaugino masses are unified at the 
GUT scale, Mqut = 3 x 10^*^ GeV and chosen a,(Mz) = 0.12. The M2 in eq.(3) is the value 
of the wino mass at the weak scale identified here with Mz- The last term in the above 
equation is a (-ve ) contribution from the D-term. It follows that the charm squark mass 
provides strong upper bound on M2 as long as m\^{MGUT) > 0: 



The branching ratio B is determined in the MSSM by the following widths |15 



T{dL ^ e+c) = ^m,, (5) 

^ ^ 167r ^ W 

(I Fl\^ + \ Fr \'){ml -ml- mlo) - Am.m^o Re{FRF*) ] (6) 

r(cL -^ xfs) = -r-^a r^^("^L'"^^'"^Y+) 

(I Gi P + I Gr HK^ - m2 - m2 + ) - 4m,m^+i?e(G^G:^) ] (7) 
where ^ Fl — '^ *^ 



2 mvK sin9\Y sinj3 

Fr = e,N, + ^ / . / iV^,, 
cosUw sinU\Y 

Q ^ msUt2 



\f2 rriw cosP' 
Gr = Vkl- 



We have adopted the same notation as in []T5[. From the expression for B in terms of 



the above decay widths, and the HERA constraint, eq.(2), one can solve for the allowed 
A'^21- The contours in the /i — M2 plane for different values of A'^21 ^^^ displayed in fig la 
{tan (3=1) and fig lb (tan/5 = 40). The horizontal lines in these figures show the upper 
bound on M2, eq.(4). We also display, the curves corresponding to two representative values 
of the chargino masses namely 45 and 85 GeV. The later is the present experimental bound 
obtained assuming R conservation. This need not hold in the presence of R violation. It 
is seen from fig. lb that for chargino mass around 85 GeV, the bound on M2 by itself rules 
out charm squark interpretation for large tan j3 independent of the value of A']^2i IIU- But 
irrespective of the value of tan/3 and the chargino mass one needs very large A'^21 — 0-13 in 
order to satisfy the bound on M2 coming from the charm squark mass. This strong bound on 
A'^21 arises because of the following reason. For M2 < 74 GeV, at least one of the charginos 
is sufficiently light and contributes dominantly to the cl decay. This reduces B [|T^ and 
results in large value for A']^2i due to eq.(2). In contrast, the chargino decay is suppressed 
kinematically for tan /3 ~ 1 if M2 > 200 GeV. This results in smaller allowed value as seen 
from the figure. But these are in conflict with the charm squark mass. 



Let us now see if one could make large A']^2i consistent with other constraints. The 



strong constraints come from atomic parity violation |T^, the decay K^ -^ tiuu [|T1|] and 
the electron neutrino mass |12[. The recent data from Cs on the relevant weak charge have 
been argued PJT8|] to imply 

A;2i < 0.074 (8) 

at 2o" level in conflict with the large value required here. In principle, the extra contribution 
due to charm squark to atomic parity violation can be canceled by a similar contribution from 
the scalar bottom or strange squark but the existing constraints on the relevant couplings 
make this cancelation difficult [|18[. Thus, one cannot easily avoid the atomic parity violation 
constraint strictly in the MSSM but this can be done by postulating new physics, e.g. the 



presence of an extra Z [18 



The other significant constraint comes from the decay K^ -^ iwv which implies [^ 

/ Tfl ~ \ 



The electron neutrino mass also gives similar constraint in the same parameter range Wl 



The question of choice of the basis becomes relevant in the discussion of these constraints. 
This is particularly so when one assumes only one A^ jr. to be non-zero. For A']^2i defined in 
the mass basis as in eq.(l) the above constraint is unavoidable if rest of the couplings are 
zero. This basis choice is natural from the point of view of interpreting HERA results but 
is not unique. One may redefine the couplings as 

and rewrite eq.(l) as follows: 

Wr = \'ijk{-Vidj + eiKl^ui)dl (10) 

HERA result would now require A'121 to be large. If this is the only non-zero X'ijk then there 
will not be any constraint on A'121 from the neutrino mass or from the K^ -^ -nvv decay 



TT| . But eq.(lO) will now generate a contribution to the neutrinoless double beta decay 



which is also severely constrained. Specifically, one has [|I3 



This clearly does not allow A'121 of 0(0. 1). Thus, notwithstanding basis dependence one has 
problem in accommodating large value for the relevant coupling. An alternative is to allow 
more than one non-zero X'ijk- It is seen from eq.(lO) that A'lji (j = l,2,3) contribute to the 
neutrinoless /3/3 decay and simultaneous presence of these may lead to cancelations. Eq.(ll) 
now gets replaced by 

(A'ln + Kl^y 121 + ^lU'm) < 2.2 x IQ-^ ( """^ V ( ""^ V^' (12) 

V 111^ 12 121^ 13 131; _ V200GeV/ V200GeVy ^ ^ 

With A'i2i ~ 0.13, cancelation between the last two terms is unlikely as it requires A'131 ~ 2. 
The first two terms can cancel but the A'm is independently constrained from the neutrino 



mass [O]. Its presence generates a large contribution to the electron neutrino mass induced 



through neutrino-gaugino mixing [O. This is given by 



o2 
m^ <i>>^ (13) 

''TT'SUSY 

The value of the induced sneutrino vev is sensitive to the MSSM parameters but can be 
approximately written as [Q 

< z/ >~ iii ma In — ^ (14) 

Requiring m^ < 2eV leads for msusY ~ 100 GeV to 

Am < .04 

It is seen that cancelations between the first two terms in eq.(12) are feasible and can allow 
A121 ~ 0.13 if this fine tuning is accepted. It must be added that the bound in the previous 
equation is quite sensitive to the MSSM parameters and for a large range in these parameters, 
the actual bound can be stronger [Q than the generic bound displayed above. 
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While wino and zino control the decay of the charm squark, its mass is mainly controlled 
by the large radiative corrections driven by the gluino mass. The unification of the gaugino 
mass parameters relates the two and leads to the above difficulty. Thus giving up this 
unification may open up a possibility of reconciling HERA events . Let us treat the gaugino 
masses Mi 2,3 at Mz as independent parameters. Then integration of the RG equation for 
the charm-squark from Mqut to Qo = 200 GeV leads to 

m^JQo) ~ ml^iMauT) + 0.77M| + 0.70M^ + OMAM^ 

+ 1/2 M| cos 2/3(1- 4/3 sin^ Ow) (15) 

If gaugino masses were to be unified at Mqut then M3 ~ 3.25M2 and Mi ~ O.5M2. Even 
in the absence of such unification, the physical gluino mass nig ~ (1 + 4.2— )M3 must be 
greater than the charm squark mass if large A121 is to be avoided. This follows since in 
the converse case, the charm squark would predominantly decay to a gluino and a quark. 
This decay being governed by strong coupling, would dominate the other decays and would 
reduce B. The later is given in case of 171^^ ^ nig by 

B ^ 1^^ ^ 2.5 X 10-3 {^\ 
32 TT a, V 0-1 / 

Such a tiny value of B would need unacceptably large A'^gi. It therefore follows that one 
must suppress the squark decay to gluino kinematically by requiring mg > tjIcj^ . Given this 
bound on M3 it follows from eq.(15) that 

M2 < 170 GeV (16) 

if nic^ ~ 220 GeV. This bound on M2 is weaker than the one in the case of the gaugino mass 
unification, eq.(4). But it nevertheless cannot suppress the decay of squarks to chargino 



kinematically. It follows |]T9| from Fig. 1 that one now approximately needs A'^gi > 0.08. 
This value is close to the 2o" limit coming from the atomic parity violation but one would 
still need some cancelations to satisfy other constraints as discussed above. Thus giving up 
unification helps only partially. 



An alternative possibility is to allow for a -ve (mass)^ for the charm squark at the 
unification scale. In view of the large positive contribution induced by the gluino mass such 
negative (mass)^ need not lead to colour breaking and may be consistent phenomeno logically. 
In fact a -ve (mass)^ for top squark has been considered in the literature |2^ in a different 



context. The universality is a simplifying feature of MSSM but it does not follow from any 
general principle. It does not hold in a large class of string based models which may allow 



negative (mass) for some sfermions as well ^^. Such masses can also arise when SUSY is 



broken by an anomalous U{1) ||2^ with some of the sparticles having -ve charge under this 
f/(l). 

The large radiative corrections induced through the running in squark masses from a high 
~ Mgut to the weak scale has played an important role in this analysis. In contrast to the 
supergravity based models, this running is over a much smaller range in models with gauge 
mediated supersymmetric breaking. But in these models, the initial value of the charm 
squark (mass)^ is positive and large with the result that these models are incompatible with 
the charm squark interpretation of HERA anomaly even without the radiative corrections 



The interpretation of HERA events in terms of stop may not suffer from the above 
mentioned difficulty encountered for the charm interpretation for two reasons. Firstly, the 
stop mass is reduced compared to the charm squark mass due to the possible large i^ — iji 
mixing as well due to the large top coupling. Secondly, this mass also involves one more 
parameter (the trilinear couplings A) compared to the charm squark mass. Thus while this 
is a less constrained possibility, imposition of the requirement that rrif: ~ 200 GeV would 
certainly lead to more constrained parameter space than considered in model independent 
studies 0. 

In summary, we have shown that the charm squark interpretation of HERA events is 
possible only for large A']^2i ~ O(O.l) in a large class of supersymmetric standard models 
characterized by a positive charm squark (mass)^ at the GUT scale. The simplest and 
the most popular minimal supergravity model with universal boundary condition falls in 
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this class. The required large value of R violating parameter is difficult to admit without 
postulating new physics and /or fine tuned cancelations. 
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Fig 1 a: The contours (continuous lines ) of constant A'121 obtained by imposing HERA 
constraint, eq.(2). The contours are for values 0.05, 0.08, 0.12, and 0.13. The horizontal 
dashed line represents the bound on M2 coming from requiring m^^ = 220 GeV. The vertical 
dash-dot lines represent the bounds on the chargino mass, the upper one for a mass of 85 
GeV and the lower one for a mass of 45 GeV. All the above are computed for tanf3 = 1. 
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Figlb 
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Fig 1 b: Same as fig la but for tanp = 40 and \[^^=0.05, 0.08, 0.12, and 0.135. 
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